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We prove ferromagnetic ordering of energy levels for XXX Heisenberg chains with spins of arbitrary 
magnitude, thus extending our previous result for the spin 1/2 chain. Ferromagnetic ordering means 
that the minimum energies in the invariant subspaces of fixed total spin are monotone decreasing 
as a function of the total spin. This result provides a ferromagnetic analogue of the well-known 
theorem by Lieb and Mattis about ordering of energy levels in antiferromagnetic and ferrimagnetic 
systems on bipartite graphs. 



PACS numbers: 75.10. Jm,75.10.Pq,75.30.Ds 

INTRODUCTION AND MAIN RESULT 

A famous theorem of Lieb and Mattis , with a subse- 
quent remark by Lieb in 0, footnote 6] , proves "ordering 
of energy levels" for a large class of Heisenberg models 
on bipartite lattices. Namely, if the two sublattices are 
A and B, and all interactions within A and B are fer- 
romagnetic while interactions in between A and B are 
antiferromagnetic, then the unique ground state multi- 
plet has total spin equal to \Sa — Sb\, where Sa and 
Sb are the maximum total spins on the two sublattices. 
Moreover, the minimum energy in the invariant subspace 
of total spin S, for S > \Sa — Sb\, is monotone increas- 
ing as a function of S. An important example is the 
usual antiferromagnet on a bipartite lattice with equal- 
size sublattices. Then the ground state is a unique spin 
singlet, and the minimum energy levels for each possible 
total spin S, are monotone increasing in 5*. 

In this paper we derive a theorem about the ordering 
of energy levels for a class of ferromagnetic Heisenberg 
models. For a ferromagnct we expect the minimum en- 
ergy for fixed total spin S to be monotone decreasing in 
S. In particular, this agrees with the well-known fact 
that the ground state has maximum total spin, and also 
with the usually tacit assumption that the first excited 
state must have total spin equal to the maximal value 
minus one. On the regular lattices this can be deduced 
by an exact spin-wave calculation. On arbitrary graphs, 
or with spins of varying magnitudes, exact results are not 
available but our argument proves that monotone order- 
ing holds at least in one-dimensional models. 

Our main object of study is the ferromagnetic Heisen- 
berg chain of L spins {L > 2). The spins may vary 
from site to site: Sx £ {5, 1, |, 2, . . .}, for a; = 1, ... , L. 



Moreover, we allow arbitrary, positive coupling constants 
Jx,x+i > 0. The quantum Hamiltonian is thus. 



H 



L-l 

^ ^ Jx.x-\-l 



SxSx-\-l 



-Sx • Sx + l 1 



(1) 



Here, S"* {i — 1,2,3) denote the standard spin-s^: matri- 
ces at site X, and Sx ■ Sx+i = J2i=i 2 3 ^x^x+i- With- 
out loss of generality, we normalize the spin matrices by 
the magnitudes, Sx, for later convenience, and subtract 
a constant so that the ground state energy vanishes. 

Since the Hamiltonian is SU(2) invariant, the vectors of 
a given total spin span an invariant subspace for H. We 
call S an admissible spin for H if this subspace contains 
at least one nonzero vector. For admissible spins, we 
define E{H, S) to be the smallest eigenvalue of H among 
all those with eigenvectors of total spin S. (We may 
define it to be infinite if S is not admissible.) Clearly, 
the largest admissible spin for H is S := J2x=i 

We say that (strict) ferromagnetic ordering of energy 
levels (FOEL) holds if it is true that E{H, S) > E{H, S') 
(with strict inequality) whenever S < S' are a pair of 
admissible spins for H. 

Theorem 1. Strict ferromagnetic ordering of energy 
levels holds for every nearest-neighbor spin chain, with 
Hamiltonian as in f7)) . 

For the special case of the spin-i chain this result was 
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proved in We conjecture that all ferromagnetic XXX 
models, i.e., on any finite graph, have the FOEL property. 
In particular we expect FOEL for the model on a lattice 
in any dimension 'l7]. In the case of constant spins, 
Sx = s, and constant couplings, Jx,x+i = J > 0, on a 
lattice, a spin-wave calculation implies that the lowest 
excitations have total spin equal to S — 1. Therefore, 
since the ground states have maximal total spin, one has 



E{H, S) > E{H, 5 - 1) > E{H, S) , 



2 



for all 5* < 5. We are not aware of any previous conjec- 
tures and certainly no proofs for other values of s or for 
the general case with arbitrary and Jx.x+i > 0. 

As a major step in the proof of Theorem 1, we will 
prove another monotonicity property of the energies 
E{H, S). Namely, suppose that H' is the Hamiltonian of 
the form Q for a ferromagnetic spin chain obtained by 
adding one or more spins to one of the ends of a spin chain 
with Hamiltonian H, and let iS' and S denote the max- 
imum total spin of these spin chains respectively. Then, 
for any < n < 5, we have 



E{H,S-n) > E{H',S' -n) 



(2) 



Note that here the deviation, n, of the maximum total 
spin is kept fixed. The inequality (0), which is part of 
Theorem 13 below, has an appealing interpretation if we 
assume that the eigenstates corresponding to the mini- 
mum energy of fixed total spin describe a droplet of op- 
posite magnetization (for the spin 1/2 XXZ chain with ei- 
ther open or periodic b.c, this has been proved in 0,0). 
(O then states that the ground state energy of such a 
droplet decreases with increasing volume, much as if the 
droplet behaves as a free particle. 



THE PROOF 

We will prove Theorem^by an induction argument in 
which the system is built up in increments of total-spin ^ . 
The total spin can be increased by i in two ways: either 
a new spin i can be appended to one end of the chain; 
or else the magnitude of the spin at the end of the chain 
can be increased by ^. We call these case (I) and case 
(II), respectively. In this way, we obtain a sequence of 
models Hi, H2, ■ ■ ■ , H^, for N = 2S, and with Lk spins 
of magnitude Sx,k, which are of the form 



Sx ' Sx-\-l 
■^x^k^x+l^k 



(3) 



H 



and such that Hi = for a single spin i, and Hj^ 
is the model of (Q. 

For any fc, if Hk+i is obtained from Hk as in case (I), 
then Hk+i = 'Hk®i 
projection Pk : Ti-k 
H'^^i = Tifc obtained from Hk, as in case (I). Then 
Pk projects onto the subspace of H'^. spanned by those 
vectors which have total spin equal to s^^ + ^ on the last 
two sites. The central part in the proof of Theorem 1 
is provided by the following theorem, which includes the 
monotonicity property mentioned in |5J. Here and later 
on, I will denote the indentity operator on C^. 

Theorem 2. Let Hi, . . . , Hf^f be a sequence of Hamil- 
tonians as described above. For each k = \, . . . , N — \, 



^ . In case (II), there is an orthogonal 
55 ^ Hk+i, and one may define 



depending on whether Hk+i is obtained from Hk as in 
case (I) or (H), one of the following two relations holds 

Hk+i>Hk®l or Hk+i = Pk (Hk ® I)P^ . (4) 

Moreover, if S is an admissible spin for Hk, then S + j 
is admissible for Hk+i, and 



E{Hk+i,S +1/2) <E{Hk,S). 



(5) 



We will prove Theorem [3 later. Assuming it, we can 
complete the proof of Theorem ^ which we will do next. 

Proof of Theorem The proof is by induction. Of 
course. Hi satisfies FOEL, simply because there is only 
one possible spin, i. Suppose that Hk satisfies strict 
FOEL. We have to show that Hk+i also satisfies FOEL. 
We will consider cases (I) and (II), separately. 

For case (I) , suppose that S < S' are such that S + ^ 
and S' + ^ are admissible for Hk+i- Note that since S' > 
S, this means that S' is admissible for Hk- Suppose that 
ip e Hk+i is the eigenvector oi Hk+i with total spin S+^, 
and energy E{Hk+i,S +5). We may assume that ■0 is a 
"highest-weight" vector, meaning that S~^ip — 0. Then, 
by the Clebsch-Gordan series for SU(2) (6(1, there are two 
highest-weight vectors f/'ijV'z G Wfe, (not both equal to 
zero) with total-spins equal to S and 5-1-1, respectively, 
and such that ip = {ipi - 2i+2^~'^'^) It) + V'2 ® ID- 
Using equation we deduce the energy inequality 

{i',Hk+ii') > (i>,{Hk®I)^P) 

= i^uHk^bi) + (1 + (25 + 2)-'){^j2,Hk^2) . 

Cross terms of the form {il^i, HkS^ 1P2) vanish since ■01 
and 02 are vectors of fixed but different total spin and 
Hk commutes with the total spin. The RHS is a con- 
vex combination of the energies of '01 and 02 with 
respect to Hk- These energies are no less than the 
minimum of E(Hk,S + 1) and E{Hk,S). But, since 
Hk satisfies FOEL, this minimum is never less than 
E{Hk,S') because S < S+l< S'. Hence, E{Hk,S) < 
E{Hk+i, S + ^). However, then using ^ we deduce that 
E{Hk+i ,S' + ^)< E{Hk+i ,S+^) which estabhshes the 
induction step. 

In case (II), we use the previous argument and the 
Rayleigh-Ritz variational principle. Namely, from I^J, we 
conclude that Hk+i — PkH'^j^^P^ , where = Hk®l 

is a Hamiltonian on Tijj^^. By the previous argument, we 
conclude that E{H'^^^,S + \) > E{Hk,S'). But by the 
variational principle, E{S + ^,Hk+i) > E{S + ^,H'i^^^) 
since the Hamiltonians are the same, but in E(S + 
^,Hk+i) one minimizes over a smaller subspace: the 
range of Pk+i. Therefore, E{Hk+i,S + i) > E{Hk,S'), 
and the rest of the argument follows case (I). This con- 
cludes the proof of Theorem 

Proof of Theorem O Q can be seen as follows. In 
case (I), where Hk+i is obtained by appending a spin-i, 
it is clear that Hk+i — Hk (23 I > since the difference 



is a Heisenberg interaction involving the new spin, which 
is non-negative by construction. In case (II), a simple 
calculation shows iJfc+i ~ P^{Hi, (g) I)Pfe. Our choice of 
spin normalization in is precisely so as to make this 
relation hold without additional numerical factors. 

To prove ©, we first observe that E{Hk,S) can be 
calculated in the sector of highest weight vectors of total 
spin S. These are vectors tp G Hk, of total spin S, such 
that S'^ip — 0, or equivalently Sf^^ip = Sip. We will cal- 
culate the matrix of Hk restricted to this subspace in a 
special basis. In this basis, we find that all off-diagonal 
matrix elements are non-positive. Let d{k, S) denote the 
dimension of the highest weight space of spin S for the 
system Hk , which is also equal to the multiplicity of the 
spin-S' representation in the chain with Hilbert space Hk ■ 
These multiplicities can be calculated by repeated use 
of the Clebsch-Gordan series for decomposing the ten- 
sor product of two spin representations into irreducible 
components. In doing this for the chains k and A: -I- 1, 
in each case starting by reducing the first two spins of 
the chain, then adding the third and so on, one sees that 
each spin-S* representation of the chain k, gives rise to a 
spin-(S'-|- 5) block in the chain k + 1. In case (I) this is so 
because the last step is adding a spin ^ . In case (II) this 
is true because the last spin is of magnitude ^ greater 
in the chain A: -I- 1 than in the chain k. Hence, we have 
d{k, S) < d{k + 1,5*+ i). Moreover, it will be possible 
to choose the two bases in such a way that 

{Hk+i\j < {Hk\j for all l<i,i< d{k, S) . 

A Perron-Frobenius style theorem 0, states that for 
any two square matrices A = (a^j), of size m, and 
B = {bki), of size n, such that: (i) n > m, (ii) all off- 
diagonal matrix elements of A and B are non-postive, 
(iii) bij < Qij for all 1 < i, j < m, one has E{B) < E{A) 
where, as before, E denotes the lowest eigenvalue. The 
minimum eigenvalues are guaranteed to be simple if there 
exist c > and an integer power p > 1 such that {cI—B)p 
has all strictly positive matrix elements. If in addition 
either to < n, or at least one of the inequalities between 
the matrix elements is strict, then E{A) > E{B). The 
inequality © is then obtained by applying this result 
with A — Hk and B = Hk+i- 

Now, we construct the basis of the highest weight 
space. One can think of the state space at site x as 
the symmetric part of 2sx spins-^. As the order of the 
spins i is irrelevant, we can label the 2s ^ + 1 states by 
the Ising configurations 

ITT---T), liT---T), liiT---T), 

E.g., liittt) is the vector normally labelled as \j,m) = 
|5/2,l/2), and nof the tensor IDfg'lb^)!!)®!?)®!?)- The 
states for a chain of L spins of magnitudes si, . . . , sl are 
then tensor products of these configurations. We shall 
call such vectors ordered Ising configurations. 
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FIG. 1: Construction of a basis vector from an ordered Ising 
configuration 

Now, we construct a basis for the spin chain consisting 
entirely of simultaneous eigenvectors of the total spin and 
its third component, with eigenvalues S and M, respec- 
tively. Start from any ordered Ising configuration such 
that 2M = # t -# j. Then, look for the leftmost j that 
has a t to its left, and draw an arc connecting this | to 
the rightmost t, left of it. At this point, one may ignore 
the paired spins, and repeat the procedure until there is 
no remaining unpaired | with an unpaired f to its left. 
This procedure guarantees that no arcs will cross and no 
arc will span an unpaired spin. The result, when ignor- 
ing all paired spins, is an ordered Ising configuration of 
a single spin. See Figure ^ for an example of this proce- 
dure. The total-spin S is equal to S minus the number of 
pairs. Clearly, the highest weight vectors are then those 
that have no unpaired J,. The vectors can be expanded 
in the tensor product basis by the following procedure: 
each arc is replaced by the spin singlet |t) fX" |i) — |i) fX" |T)i 
and the unpaired spins are replaced by their tensor prod- 
ucts. Finally, one symmetrizes in each block. In the case 
Sx ~ 1/2 for all x, this was called the Hulthcn bracket 
basis by Temperley and Lieb |^ , who also extended it to 
SUg(2) and proved that it is a basis. Frenkel and Kho- 
vanov introduced the general case 0. 

Next, we calculate the matrix representation of the 
Hamiltonians Hk in the invariant subspaces of all highest 
weight vectors of a given total spin. This is most easily 
accomplished by deriving a graphical representation for 
the action of each term in the Hamiltonian. We start with 
the case where all Sx — 1/2. Then the interaction terms 
generate the Temperley-Lieb algebra with 9 = 1 d, ^3 ■ 
The action of the negative of the interaction term in the 
Hamiltonian can be graphically represented by 

1 0. P 

— hx,x+l — 7: (45a: • Sx + 1 — 1) = 

2 □ Q 

We refer to the graphical diagram for —hx^x+i as Ux.x+i- 
Let a be the diagram for the basis vector \a). The graph- 
ical rules are as follows: (i) if x and x + 1 are both un- 
paired in a, we have hx,x+i\c^) — 0, (ii) if the composition 
of a and Ux,x+i is isotopic to (3, with (3 ^ a, then the 
hx,x+i\a) = -|/3), (in) if a = /3, then hx,x+i\a) = 2\a). 
This means that every loop contributes a factor of 2 to 
hx^x+i- But this only happens when the "cup" of Ux,x+i 
is paired with an arc in a. This is clearly a diagonal 
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term. If the last site of an arc for a is x, then there is no 
loop for hy^y+i with y > x. 

To generalize this to arbitrary values of Sx, it is suf- 
ficient to write the Heisenberg interaction for arbitrary 
spins as an interaction between spin i 's making up the 
spin Sx and Sx+i, conjugated with the projections onto 
the symmetric vectors. 

The result is the following: 



-hx,x+i — 7:( Sx ■ Sx+i — 1) — 

2 SxSx+l 



2sx 2sx+i 



2sx 2sx+i 



Here, the rectangles with label 2s represent the sym- 
metrizing projections on the space of 2s spin ^ variables. 
The fundamental algebraic property that will allow us 
to calculate the matrix elements of Hk graphically is the 
Jones- Wenzl relation (c.f., IflJ and references therein): 



\2s+l 




2s 



2s + 1 



(6) 



For any element of the basis introduced above one can 
compute the action of the Hamiltonian and write it as a 
linear combination of the same basis vectors. From the 
grahical rules it is easy to observe that all off-diagonal 
matrix elements are non-positive. 

Furthermore, with each basis vector for the system Hk, 
we can identify a basis vector for the system Hk+i- First, 
we consider the case where one spin 1/2 has been added 
to the system Hk- In that case one can simply add one 
box with one up-spin to each basis vector. This then cor- 
responds to a basis vectors for the system where 
there is no arc to the last spin. The case of increasing 
the magnitude of one of the spins by 1 /2 can be treated 
in a similar way by again using the Jones- Wenzl projec- 
tion. In this case, the label of the rightmost box in any 
basis vector for the Hk system is raised by one but the 
number of arcs remains unchanged. Again, this clearly 
leads to a subset of the basis vectors for the Hk+i sys- 
tem. The crucial property that allows us to compare the 
two Hamiltonians is the following. When H^+i acts on 
a basis vector obtained from a corresponding H^ vector 
as we have just described, the only possible new terms 
that are generated are off-diagonal terms, which do not 
contain a bubble and, hence, are negative. The details 
of the calculation of these matrix elements and further 
applications will appear elsewhere 



EXTENSIONS AND DISCUSSION 

Our results can be extended in various directions. We 
discuss two generalizations in particular. 



The first is the analogous problem for the anisotropic 
model, i.e., the ferromagnetic XXZ model, with 
anisotropy A > 1. Of course, the XXZ Hamiltonian no 
longer has the full SU(2) symmetry, but for the spin-^ 
chain there is an SUq(2)-symmetric model obtained just 
by adding suitable boundary fields 0|. The represen- 
tation theory of SUg(2) is similar to that of SU(2)|0|. 
Using the SUg(2) symmetry, one can show that, for a 
chain of L sites, the first excited state has total "spin" 
[L — l)/2 By replacing the Hulthen brackets with 

their g— analogues, FOEL can be proved for the spin-i 
XXZ chains j^]- In principle one can generalize this re- 
sult to higher spin models as well, but one has to replace 
the standard XXZ interaction by a quantum group sym- 
metric interaction. Since such examples are of a lesser 
physical relevance we do not pursue this. 

The second type of generalization concerns higher or- 
der interactions with the full SU(2) symmetry. As our 
argument to prove Theorem ^ only uses the symmetry of 
the Hamiltonian, and not its exact form, that result car- 
ries over immediately. The dual canonical basis of highest 
weight vectors also does not depend on the Hamiltonian. 
Theorem |5] requires, however, that the off-diagonal ele- 
ments of the Hamiltonian with respect to this basis are 
non-negative. This puts a restriction on the class of in- 
teractions. As an example, we consider the spin-1 chain. 
For this case there is one non-trivial parameter in the 
family of SU(2)-invariant nearest neighbor interactions, 
say the coefficient of the biquadratic term: 



H = Y,{l-Sx 



Sx+l) + ^(1 — {Sx ■ Sx + l) 



(7) 



It can be proved that the Hamiltonian lO has FOEL 
whenever < i < 1/3, and this result is optimal. 

As a practical application of our results we point out 
the following. Let H he & Hamiltonian satisfying FOEL 
and suppose one would like to find all the eigenvalues of 
H below energy E. It is sufficient to diagonalize H in 
the subspaces with total spin S ~ k, for k = 0,1, .. .K, 
such that E{H, S-K) <E. The first value K such that 
E{H, S — K — 1) > E \s large enough. It is then guaran- 
teed that one has all eigenvalues of the full H below that 
value. This is interesting because the dimension of the 
subspace for a chain of length L with total spin equal to 
S — k\s 0{L^), while the full Hilbert space has dimension 
(25' -I- 1)^ (when all particles have spin S). In particular, 
one can obtain the gap by considering K — 1. 
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